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Detrended fluctuation analysis (DFA) is a scaling analysis method used to estimate long-range 
power-law correlation exponents in noisy signals. Many noisy signals in real systems display trends, 
so that the scaling results obtained from the DFA method become difficult to analyze. We system- 
atically study the effects of three types of trends — linear, periodic, and power-law trends, and offer 
examples where these trends are likely to occur in real data. We compare the difference between the 
scaling results for artificially generated correlated noise and correlated noise with a trend, and study 
how trends lead to the appearance of crossovers in the scaling behavior. We find that crossovers 
result from the competition between the scaling of the noise and the "apparent" scaling of the trend. 
We study how the characteristics of these crossovers depend on (i) the slope of the linear trend; 
(ii) the amplitude and period of the periodic trend; (iii) the amplitude and power of the power-law 
trend and (iv) the length as well as the correlation properties of the noise. Surprisingly, we find that 
the crossovers in the scaling of noisy signals with trends also follow scaling laws — i.e. long-range 
power-law dependence of the position of the crossover on the parameters of the trends. We show 
that the DFA result of noise with a trend can be exactly determined by the superposition of the 
separate results of the DFA on the noise and on the trend, assuming that the noise and the trend 
are not correlated. If this superposition rule is not followed, this is an indication that the noise 
and the superimposed trend are not independent, so that removing the trend could lead to changes 
in the correlation properties of the noise. In addition, we show how to use DFA appropriately to 
minimize the effects of trends, and how to recognize if a crossover indicates indeed a transition from 
one type to a different type of underlying correlation, or the crossover is due to a trend without any 
transition in the dynamical properties of the noise. 



I. INTRODUCTION 

Many physical and biological systems exhibit com- 
plex behavior characterized by long-range power-law cor- 
relations. Traditional approaches such as the power- 
spectrum and correlation analysis are not suited to accu- 
rately quantify long-range correlations in non-stationary 
signals — e.g. signals exhibiting fluctuations along poly- 
nomial trends. Detrended fluctuation analysis (DFA) 
[Q-Q is a scaling analysis method providing a simple 
quantitative parameter — the scaling exponent a — to 
represent the correlation properties of a signal. The ad- 
vantages of DFA over many methods are that it per- 
mits the detection of long-range correlations embedded 
in seemingly non-stationary time series, and also avoids 
the spurious detection of apparent long-range correla- 
tions that are artifact of non-stationarity. In the past 
few years, more than 100 publications have utilized the 
DFA as method of correlation analysis, and have uncov- 
ered long-range power-law correlations in many research 
fields such as cardiac dynamics |fi|~p3l|, bioinformatics 
^|,3, |2^ |3j3l , economic s [j35r - f47| ], meteorology |4g| |50|], ge- 
ology []5 lfi~7 ethology [|52fl etc. Furthermore, the DFA 
method may help identify different states of the same 
system according to its different scaling behaviors — e.g. 
the scaling exponent a for heart inter-beat intervals is 
different for healthy and sick individuals (14,10, 17[p3|. 

The correct interpretation of the scaling results ob- 
tained by the DFA method is crucial for understanding 



the intrinsic dynamics of the systems under study. In 
fact, for all systems where the DFA method was applied, 
there are many issues that remain unexplained. One 
of the common challenges is that the correlation expo- 
nent is not always a constant (independent of scale) and 
crossovers often exist — i.e. change of the scaling expo- 
nent a for different range of scales §,||,||. A crossover 
usually can arise from a change in the correlation proper- 
ties of the signal at different time or space scales, or can 
often arise from trends in the data. In this paper, we sys- 
tematically study how different types of trends affect the 
apparent scaling behavior of long-range correlated sig- 
nals. The existence of trends in times series generated by 
physical or biological systems is so common that it is al- 
most unavoidable. For example, the number of particles 
emitted by a radiation source in an unit time has a trend 
of decreasing because the source becomes weaker p4],p5[; 
the density of air due to gravity has a trend at different 
altitude fl56fl ; the air temperature in different geographic 
locations and the water flow of rivers have a periodic 
trend due to seasonal changes p9|,[50|,p7|-[59|]; the occur- 
rence rate of earthquakes in certain area has trend in 
different time period p0| |. An immediate problem fac- 
ing researchers applying scaling analysis to time series 
is whether trends in data arise from external conditions, 
having little to do with the intrinsic dynamics of the sys- 
tem generating noisy fluctuating data. In this case, a 
possible approach is to first recognize and filter out the 
trends before we attempt to quantify correlations in the 
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noise. Alternatively, trends may arise from the intrinsic 
dynamics of the system, rather than being an epiphe- 
nomenon of external conditions, and thus may be corre- 
lated with the noisy fluctuations generated by the sys- 
tem. In this case, careful considerations should be given 
if trends should be filtered out when estimating correla- 
tions in the noise, since such "intrinsic" trends may be 
related to the local properties of the noisy fluctuations. 

Here we study the origin and the properties of 
crossovers in the scaling behavior of noisy signals, by ap- 
plying the DFA method first on correlated noise and then 
on noise with trends, and comparing the difference in the 
scaling results. To this end, we generate artificial time 
series — anticorrelated, white and correlated noise with 
standard deviation equal to one — using the modified 
Fourier filtering method introduced by Makse et al. |33) . 
We consider the case when the trend is independent of 
the local properties of the noise (external trend). We find 
that the scaling behavior of noise with a trend is a su- 
perposition of the scaling of the noise and the apparent 
scaling of the trend, and we derive analytical relations 
based on the DFA, which we call "superposition rule" . 
We show how this "superposition rule" can be used to 
determine if the trends are independent of the noisy fluc- 
tuation in real data, and if filtering these trends out will 
no affect the scaling properties of the data. 

The outline of this paper is as follows. In Sec.|l[ we re- 
view the algorithm of the DFA method, and in Appendix 
[a| we compare the performance of the DFA with the clas- 
sical scaling analysis — Hurst's analysis (R/S analysis) — 
and show that the DFA is a superior method to quan- 
tify the scaling behavior of noisy signals. In Sec. [II, 



we consider the effect of a linear trend and we present 
an analytic derivation of the apparent scaling behavior 
of a linear trend in Appendix JQ. In Sec. IV, we study 
a periodic trend, and in Sec. fVTthe effect of power-law 
trend. We systematically study all resulting crossovers, 
their conditions of existence and their typical character- 
istics associated with the different types of trends. In 
addition, we also show how to use DFA appropriately to 
minimize or even eliminate the effects of those trends in 
cases that trends are not choices of the study, that is, 
trends do not reflect the dynamics of the system but are 
caused by some "irrelevant" background. Finally, Sec. VI 
contains a summary. 
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and is divided into boxes of equal size, n. In each box, we 
fit the integrated time series by using a polynomial func- 
tion, yfit(i), which is called the local trend. For order-£ 
DFA (DFA-1 if £ = 1, DFA- 2 if £ = 2 etc.), I order poly- 
nomial function should be applied for the fitting. We 
detrend The integrated time series, y(i) by subtracting 
the local trend yfit(i) in each box, and we calculate the 
detrended fluctuation function 



=v{i)-Vfit(i). 
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For a given box size n, we calculate the root mean square 
(rms) fluctuation 
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The above computation is repeated for box sizes n (dif- 
ferent scales) to provide a relationship between F(n) and 
n. A power-law relation between F(n) and the box size 
n indicates the presence of scaling: F(n) ~ n a . The 
parameter a, called the scaling exponent or correlation 
exponent, represents the correlation properties of the sig- 
nal: if a = 0.5, there is no correlation and the signal is an 
uncorrelated signal (white noise); if a < 0.5, the signal is 
anticorrelated; if a > 0.5, there are positive correlations 
in the signal. 



III. NOISE WITH LINEAR TRENDS 

First we consider the simplest case: correlated noise 
with a linear trend. A linear trend 



u(i) = A L i 



(5) 



is characterized by only one variable — the slope of the 
trend, A^. For convenience, we denote the rms fluctu- 
ation function for noise without trends by F v (n), linear 
trends by FL in) , and noise with a linear trend by F^l (n) ■ 



II. DFA 

To illustrate the DFA method, we consider a noisy time 
series, u(i) (i = 1, .., N max ). We integrate the time series 
u(i), 
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A. DFA-1 on noise with a linear trend 

Using the algorithm of Makse (63|, we generate cor- 
related noise with standard deviation one, with a given 
correlation property characterized by a given scaling ex- 
ponent a. We apply DFA-1 to quantify the correlation 
properties of the noise and find that only in certain good 
fit region the rms fluctuation function F„(n) can be ap- 
proximated by a power-law function [see Appendix |a| 

F v (n) = b n a (6) 
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where bo is a parameter independent of the scale n. We 
find that the good fit region depends on the correlation 
exponent a [see Appendix |X| . We also derive analyti- 
cally the rms fluctuation function for linear trend only 
for DFA-1 and find that [see Appendix ^] 



F L (n) = k A h n c 



(7) 



where kg is a constant independent of the length of trend 
N max , of the box size n and of the slope of the trend A^. 
We obtain otL = 2. 
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FIG. 1. Crossover behavior of the root mean square fluc- 



tuation function F^fa) for noise (of length N„ 



and 



correlation exponent a — 0.1) with superposed linear trends 
of slope Ah — 2~ 16 , 2~ 12 , 2~ 8 . For comparison, we show F v (n) 
for the noise (thick solid line) and Fi,[n) for the linear trends 
(dot-dashed line) (Eq.dW)). The results show that a crossover 
at a scale n x for F n ^(n). For n < n x , the noise dominates 
and F v i,{n) ~ F v (n). For n > n x , the linear trend domi- 
nates and Frfhiji) ~ Fhin). Note that the crossover scale n x 
increases when the slope Al of the trend decreases. 

Next we apply the DFA-1 method to the superposi- 
tion of a linear trend with correlated noise and we com- 
pare the rms fluctuation function F^in) with F^in) [see 
Figjl[. We observe a crossover in F^n) at scale n = n x . 
For n < n x , the behavior of F V L(n) is very close to the 
behavior of Fi,(n), while for n > n x , the behavior of 
-Ft)l(^) is very close to the behavior of Fl(ti). A sim- 
ilar crossover behavior is also observed in the scaling 
of the well-studied biased random walk ]6l],|6^|. It is 
known that the crossover in the biased random walk is 
due to the competition of the unbiased random walk and 
the bias [see Fig. 5. 3 of [p2|]. We illustrate this observa- 
tion in Fig. ^, where the detrended fluctuation functions 
(Eq. (||)) of the correlated noise, Y v (i), and of the noise 
with a linear trend, Y^li) are shown. For the box size 
n < n x as shown in Fig. ||(a) and (b), Y^li) « Y v (i). 
For n> n x as shown in Fig. ||(c) and (d), Y r) -^(i) has dis- 



tinguishable quadratic background significantly different 
from Y v (i). This quadratic background is due to the inte- 
gration of the linear trend within the DFA procedure and 
represents the detrended fluctuation function Yl of the 
linear trend. These relations between the detrended fluc- 
tuation functions Y(i) at different time scales n explain 
the crossover in the scaling behavior of F^ln): from very 
close to F n (n) to very close to F^iri) (observed in FigJ^). 
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FIG. 2. Comparison of the detrended fluctuation function 
for noise Y v (i) and for noise with linear trend Y v i^(i) at differ- 
ent scales, (a) and (c) are Y v for noise with a = 0.1; (b) and 
(d) are Y v h for the same noise with a linear trend with slope 
j4l = 2~ 12 (the crossover scale n x — 320 see Fig. ^). (a) (b) 
for scales n < n x the effect of the trend is not pronounced 
and Y v ~ Y v l (i.e. Y v 3> (c)(d) for scales n > n x , the 
linear trend is dominant and Y v <C I^l- 

The experimental results presented in Figs.|l| and ^ sug- 
gest that the rms fluctuation function for a signal which 
is a superposition of a correlated noise and a linear trend 
can be expressed as: 



[F vL (n)] 2 = [F L (n)} 2 



[F v (n)} 2 



(8) 



We provide an analytic derivation of this relation in Ap- 
pendix ||, where we show that Eq. (||) holds for the super- 
position of any two independent signals — in this particu- 
lar case noise and a linear trend. We call this relation the 
"superposition rule" . This rule helps us understand how 
the competition between the contribution of the noise 
and the trend to the rms fluctuation function F^^n) at 



different scales n leads to appearance of crossovers |61|. 

Next, we ask how the crossover scale n x depends on: 
(i) the slope of the linear trend v4l, (ii) the scaling ex- 
ponent a of the noise, and (iii) the length of the signal 
N m ax- Surprisingly, we find that for noise with any given 
correlation exponent a the crossover scale n x itself fol- 
lows a power-law scaling relation over several decades: 
n x ~ (Al) 6 (see Fig. ||). We find that in this scaling re- 
lation, the crossover exponent 9 is negative and its value 
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depends on the correlation exponent a of the noise — the 
magnitude of 9 decreases when a increases. We present 
the values of the "crossover exponent" for different cor- 
relation exponents a in Table 
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FIG. 3. The crossover n x of F^l(i-) for noise with a lin- 
ear trend. We determine the crossover scale n x based on the 
difference A between log F v (noise) and log F v l (noise with 
a linear trend). The scale for which A = 0.05 is the esti- 
mated crossover scale n x . For any given correlation exponent 
a of the noise, the crossover scale n x exhibits a long-range 
power-law behavior n x ~ (A^) 6 , where the crossover expo- 
nent 9 is a function of a [see Eq.(H) and Table n|. 



TABLE I. The crossover exponent 6 from the 
power- law relation between the crossover scale n x and 
the slope of the linear trend Al — n x ~ (Ah) 6 — for dif- 
ferent values of the correlation exponents a of the noise 
[Fig. |^ . The values of 6 obtained from our simulations 
are in good agreement with the analytical prediction 
-1/(2 - a) [Eq. (|)]. Note that -1/(2 - a) are not 
always exactly equal to 6 because F v (n) in simulations 
is not a perfect simple power-law function and the way 
we determine numerically n x is just approximated. 
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To understand how the crossover scale depends on the 
correlation exponent a of the noise we employ the super- 
position rule [Eq.(||)] and estimate n x as the intercept 
between F v (n) and F L (n). From the Eqs. @ and M), we 
obtain the following dependence of n x on a: 
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This analytical calculation for the crossover exponent 
— 1/(c*l — ol) is in a good agreement with the observed 
values of 9 obtained from our simulations [see Fig.|] and 
Table | . 

Finally, since the F^(n) does not depend on N max as 
we show in Eq.(Q) and in Appendix |^, we find that n x 
does not depend on N max . This is a special case for 
linear trends and does not always hold for higher order 
polynomial trends [see Appendix 



B. DFA-2 on noise with a linear trend 
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FIG. 4. Comparison of the rms fluctuation function F v (n) 
for noise with different types of correlations (lines) and F v l (n) 
for the same noise with a linear trend of slope Al = 2~ 12 
(symbols) for DFA-2. F v i,(n) = F v (n) because the inte- 
grated linear trend can be perfectly filtered out in DFA-2, 
thus Yh(i) = from Eq.(||). We note, that to estimate accu- 
rately the correlation exponents one has to choose an optimal 
range of scales n, where F(n) is fitted. For details see Ap- 
pendix [a| 



Application of the DFA-2 method to noisy signals with- 
out any polynomial trends leads to scaling results identi- 
cal to the scaling obtained from the DFA-1 method, with 
the exception of some vertical shift to lower values for the 
rms fluctuation function F v (n) [see Appendix ^] . How- 
ever, for signals which are a superposition of correlated 
noise and a linear trend, in contrast to the DFA-1 results 
presented in Fig. l], F^^n) obtained from DFA exhibits 
no crossovers, and is exactly equal to the rms fluctuation 
function F v (n) obtained from DFA-2 for correlated noise 
without trend (see Fig. ||). These results indicate that 
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a linear trend has no effect on the scaling obtained from 
DFA-2. The reason for this is that by design the DFA-2 
method filters out linear trends, i.e. Yl(i) = (Eq.( ||)) 
and thus F^^n) — F^fn) due to the superposition rule 
(Eq. (ph). For the same reason, polynomial trends of or- 
der lower than I superimposed on correlated noise will 
have no effect on the scaling properties of the noise when 
DFA-£ is applied. Therefore, our results confirm that the 
DFA method is a reliable tool to accurately quantify cor- 
relations in noisy signals embedded in polynomial trends. 
Moreover, the reported scaling and crossover features of 
F(n) can be used to determine the order of polynomial 
trends present in the data. 



IV. NOISE WITH SINUSOIDAL TREND 

In this section, we study the effect of sinusoidal trends 
on the scaling properties of noisy signals. For a signal 
which is a superposition of correlated noise and sinu- 
soidal trend, we find that based on the superposition rule 
(Appendix |b|) the DFA rms fluctuation function can be 
expressed as 



[F vS (n)f = [F v (n)] 2 



[F s (r 



(10) 



where F^s (n) is the rms fluctuation function of noise with 
a sinusoidal trend, and Fs(n) is for the sinusoidal trend. 
First we consider the application of DFA-1 to a sinu- 
soidal trend. Next we study the scaling behavior and the 
features of crossovers in F v s(n) for the superposition of 
correlated noise and sinusoidal trend employing the su- 
perposition rule [Eq.(|n])]. At the end of this section, we 
discuss the results obtained from higher order DFA. 



A. DFA-1 on sinusoidal trend 

Given a sinusoidal trend u(i) — As sin (2m/T) (i = 
1, N max ), where As is the amplitude of the signal and 
T is the period, we find that the rms fluctuation func- 
tion Fs(n) does not depend on the length of the signal 
N max , and has the same shape for different amplitudes 
and different periods [Fig. |) . We find a crossover at scale 
corresponding to the period of the sinusoidal trend 



n 2 x 



T, 
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and does not depend on the amplitude As- We call this 
crossover n 2x for convenience, as we will see later. For 
ri < n 2 x , the rms fluctuation Fg (n) exhibits an ap- 
parent scaling with the same exponent as Fi,(n) for the 
linear trend [see Eq. (m)]: 



F S (n) 



(12) 



where k\ is a constant independent of the length N max , 
of the period T and the amplitude As of the sinusoidal 



signal, and of the box size n. As for the linear trend 
[Eq. (0)] , we obtain as — 2 because at small scales (box 
size n) the sinusoidal function is dominated by a linear 
term. For n > n 2x , due to the periodic property of the 
sinusoidal trend, Fs(n) is a constant independent of the 
scale n : 



F S (n) 
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2V2vr 



As • T. 



(13) 



The period T and the amplitude As also affects the ver- 
tical shift of Fs(n) in both regions. We note that in 
Eqs.(p"2|) and (jl^), Fs(n) is proportional to the ampli- 
tude Ag , a behavior which is also observed for the linear 
trend [Eq. (|)]. 
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FIG. 5. Root mean square fluctuation function Fs(n) for 
sinusoidal functions of length N max — 2 17 with different am- 
plitude As and period T. All curves exhibit a crossover at 
7i2x ~ T/2, with a slope as = 2 for n < n2x, and a flat 
region for n > ri2x • There are some spurious singularities at 
n — j-j (j is a positive integer) shown by the spikes. 



B. DFA-1 on noise with sinusoidal trend 

In this section, we study how the sinusoidal trend af- 
fects the scaling behavior of noise with different type of 
correlations. We apply the DFA-1 method to a signal 
which is a superposition of correlated noise with a sinu- 
soidal trend. We observe that there are typically three 
crossovers in the rms fluctuation i 7 i 7 s(7i) at characteristic 
scales denoted by tii x , 7i2 X and 773 x [Fig. [|. These three 
crossovers divide i^g(n) into four regions, as shown in 
Fig. ^(a) (the third crossover cannot be seen in Fig. ^(b) 
because its scale n 3x is greater than the length of the sig- 
nal) . We find that the first and third crossovers at scales 
Tiix and 773X respectively [see Fig. ^| result from the com- 
petition between the effects on F„s (n) of the sinusoidal 
signal and the correlated noise. For n < n\ x (region I) 
and n > 7i3x (region IV), we find that the noise has the 
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dominating effect (F r) (n) > Fs(n)), so the behavior of 
Fris(n) is very close to the behavior of F v (n) [Eq. 
For ni x < n < ri2x (region II) and ri2x < n < «3x (re- 
gion III) the sinusoidal trend dominates (Fs(n) > F v (n)), 
thus the behavior of F v s (n) is close to Fg (n) [see Fig. ^ 
and Fig. § . 
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FIG. 6. Crossover behavior of the root mean square fluctu- 
ation function F v s(n) (circles) for correlated noise (of length 
N m ax = 2 ) with a superposed sinusoidal function charac- 
terized by period T = 128 and amplitude As = 2. The rms 
fluctuation function F v (n) for noise (thick line) and Fg (n) for 
the sinusoidal trend (thin line) are shown for comparison, (a) 
F v s(n) for correlated noise with a — 0.9. (b) F v s(n) for an- 
ticorrelated noise with a — 0.9. There are three crossovers in 
F,s(fi)i at scales nix, "2x and ri3 X (the third crossover can 
not be seen in (b) because it occurs at scale larger than the 
length of the signal). For n < nix and n > mx, the noise 
dominates and F v s(n) ~ F v (n) while for nix < n < ns x , 
the sinusoidal trend dominates and F v s(n) « Fs(n). The 
crossovers at nix and ri3 X are due to the competition between 
the correlated noise and the sinusoidal trend [see Fig. [?]] , while 
the crossover at ri2 X relates only to the period T of the sinu- 
soidal [Eq. (El)]. 
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FIG. 7. Comparison of the detrended fluctuation function 
for noise, Y v (i) and noise with sinusoidal trend, Y v s(i) in four 
regions as shown in Fig. |^. The same signals as in Fig. ^| are 
used. Panels (a)-(f) correspond to Fig. ||(b) for anticorrelated 
noise with exponent a — 0.1, and panels (g)-(h) correspond 
to the Fig. |^(a) for correlated noise with exponent a — 0.9. 
(a)-(b) For all scales n < nix, the effect of the trend is not 
pronounced and Y v s(i) ~ Y v (i) leading to F v s(n) « F v (n) 
(Fig. ^(a)). (c)(d) For n2 X > n > ni x , the trend is domi- 
nant, Kjs(i) » Y v (i) and F v s{n) « Fs(n). Since n 2x ~ T/2 
(Eq. (|l l|) ) , the scale n < T/2 and the sinusoidal behavior 
can be approximated as a linear trend. This explains the 
quadratic background in Y v s(i) (d) [see Fig. ^(c)(d)]. (e)(f) 
For 7i2x < n < 7i3x (i.e. n S> T/2), the sinusoidal trend 
again dominates — Y v s (i) is periodic function with period T. 
(g)(h) for n > n3x, the effect of the noise is dominant and 
the scaling of F v s follows the scaling of F v (Fig. pi a)). 



To better understand why there are different regions in 
the behavior of F v s(n), we consider the detrended fluc- 
tuation function [Eq. ([|) and Appendix |B| of the corre- 
lated noise YjJi), and of the noise with sinusoidal trend 
Y 7j s- In Fig. [7| we compare Y v (i) and ^s(*) for anticor- 
related and correlated noise in the four different regions. 
For very small scales n < riix, the effect of the sinu- 
soidal trend is not pronounced, Y^sty) Y^i), indicat- 
ing that in this scale region the signal can be considered 
as noise fluctuating around a constant trend which is fil- 
tered out by the DFA-1 procedure [Fig. 0(a) (b)]. Note, 
that the behavior of Y v $ [Fig. 0(b)] is identical to the be- 
havior of I^l [Fig. ||(b)], since both a sinusoidal with 
a large period T and a linear trend with small slope 
Al can be well approximated by a constant trend for 



G 



n < n% x - For small scales ni x < n < «2x (region II), we 
find that there is a dominant quadratic background for 
Y v s(i) [Fig. 0(d)]. This quadratic background is due to 
the integration procedure in DFA-1, and is represented 
by the detrended fluctuation function of the sinusoidal 
trend Yg(i). It is similar to the quadratic background 
observed for linear trend Y v ^(i) [Fig. ||(d)] — i.e. for 
nix < n < ri2x the sinusoidal trend behaves as a linear 
trend and Ys(i) w Y^(i). Thus in region II the "lin- 
ear trend" effect of the sinusoidal is dominant, Y$ > Y v , 
which leads to F n s(n) ~ Fs(n). This explains also why 
-Ft)s(^) for n < n 2x (Fig. g) exhibits crossover behav- 
ior similar to the one of F^{n) observed for noise with 
a linear trend. For ri2x < n < n^x (region III) the 
sinusoidal behavior is strongly pronounced [Fig. 0(f)], 
Ys(i) ^> Y v (i), and Y^i) w Y$(i) changes periodically 
with period equal to the period of the sinusoidal trend 
T. Since Y^ii) is bounded between a minimum and a 
maximum value, F v s(n) cannot increase and exhibits a 
flat region (Fig. ^). At very large scales, n > n 3x , the 
noise effect is again dominant (Ys(i) remains bounded, 
while Y n grows when increasing the scale) which leads to 
i^s(n) « Frjin), and a scaling behavior corresponding to 
the scaling of the correlated noise. 
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(a) Noise + sin. trend (A s =5.0) 




(b) Noise + sin. trend (T=2 ) 
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(d) Noise + sin. trend (A s =2) 
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(e) Noise + sin. trend (T=16) 
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FIG. 8. Dependence of the three crossovers in F v s(n) for 
noise with a sinusoidal trend (Fig. |^) on the period T, and 
amplitude As of the sinusoidal trend, (a) Power-law rela- 
tion between the first crossover scale ni x and the period T 
for fixed amplitude As and varying correlation exponent a: 
riix ~ T 0T1 , where Oti is a positive crossover exponent [see 
Table |fl] and Eq. [lij. (b) Power- law relation between the 
first crossover nix an d the amplitude of the sinusoidal trend 
As for fixed period T and varying correlation exponent a: 
nix ~ Ag A1 where Oai is a negative crossover exponent [Ta- 
ble [n| and Eq. (0)]. (c) The second crossover scale ii2x de- 
pends only on the period T: nix ~ j iS t2 j wnere q T2 ~ i. 

(d) Power-law relation between the third crossover 713 x and 
T for fixed amplitude As and varying a trend: 713 x ~ T 9 " 13 . 

(e) Power-law relation between the third crossover ri3 X and 

^ A3 . We find that 



ll',y 



As for fixed T and varying a 
8as = 0T3 [Table |[li| and Eq. © 



(^; 



First, we consider ni x . Surprisingly, we find that for 
noise with any given correlation exponent a the crossover 
scale nix exhibits long-range power-law dependence of 
the period T — ni x ~ T 0T1 , and the amplitude As - 
nix ~ (As) 6ai of the sinusoidal trend [see Fig. ^(a) and 
(b)]. We find that the "crossover exponents" #xi and 
#ai have the same magnitude but different sign — 9ti is 
positive while Oai is negative. We also find that the mag- 
nitude of 8ti and 9ai increases for the larger values of 
the correlation exponents a of the noise. We present the 
values of 6>xi and 8ai for different correlation exponent 
a in Table |n[ To understand these power-law relations 
between n\ x and T, and between ni x and As, and also 
how the crossover scale ni x depends on the correlation 
exponent a we employ the superposition rule [Eq. [10] 
and estimate ni x analytically as the first intercept n^x 
of F v (n) and Fs(n). From Eqs. ( |l2| ) and (g), we obtain 
the following dependence of ni x on T, As and a: 



hi. 

ki A s 



l/(2-a) 



(14) 



From this analytical calculation we obtain the fol- 
lowing relation between the two crossover expo- 
nents #ti and 9ai and the correlation exponent a: 
T1 = — Q A1 = 1/(2 — cv), which is in a good agree- 
ment with the observed values of £?xi , #ai obtained from 
simulations [see Fig. |(a) (b) and Table ||. 

Next, we consider n2 X . Our analysis of the rms fluc- 
tuation function Fs (n) for the sinusoidal signal in Fig. || 
suggests that the crossover scale Fs(n) does not depend 
on the amplitude As of the sinusoidal. The behavior of 
the rms fluctuation function F n s (n) for noise with super- 
imposed sinusoidal trend in Fig. ||(a) an d (b) indicates 
that n2x does not depend on the correlation exponent 
a of the noise, since for both correlated (a = 0.9) and 
anticorrelated (a = 0) noise (T and As are fixed), the 
crossover scale n2 X remains unchanged. We find that n2x 
depends only on the period T of the sinusoidal trend and 
exhibits a long-range power-law behavior n2 X ~ 
with a crossover exponent 0X2 ~ 1 (Fig. ||(c)) which is in 



agreement with the prediction of Eq. ( |1 1| ) . 

For the third crossover scale n3 X , as for ni x we find 
a power-law dependence on the period T, n3 X ~ t 9t3 , 
and amplitude As, n^x ~ (As) A3 , of the sinusoidal trend 
[see Fig. ||(d) and (e)]. However, in contrast to the m x 
case, we find that the crossover exponents 9t P 3 and 6*A3 
are equal and positive with decreasin g va lues for increas- 
ing correlation exponents a. In Table III, we present the 



values of these two exponents for different correlation ex- 
ponent a. To understand how the scale n3 X depends on 
T, As and the correlation exponent a simultaneously, 
we again employ the superposition rule [Eq. (|To|)] and 
estimate n3 X as the second intercept nMl of F„(n) and 



Fs(n). From Eqs. 
dependence: 



(|13|) and (o), we obtain the following 



n 3x = 



\2V2nbo 



AsT 



l/a 



(15) 



From this analytical calculation we obtain 0T3 = $A3 = 
l/a which is in good agreement with the values of 8t3 
and 9 A3 observed from simulations [Table III ] . 



TABLE II. The crossover exponents Oti and Oai 
characterizing the power- law dependence of nix on the 
period T and amplitude As obtained from simulations: 
nix ~ T 0T1 and nix ~ (As) 9ai for different value 
of the correlation exponent a of noise [Fig. ^(a)(b)]. 
The values of Oti and 9ai are in good agreement with 
the analytical predictions 8ti = — Oai = 1/(2 — a) 
[Eq. @]. 



a 


Oti 


-0A1 


1/(2 -a) 


0.1 


0.55 


0.54 


0.53 


0.3 


0.58 


0.59 


0.59 


0.5 


0.66 


0.66 


0.67 


0.7 


0.74 


0.75 


0.77 


0.9 


0.87 


0.90 


0.91 



TABLE III. The crossover exponents 0t3 and 0a3 
for the power- law relations: 713 x ~ y fl T3 anc j 
n3x ~ (As)° A3 for different value of the correlation 
exponent a of noise [Fig. |^(c)(d)]. The values of P 3 
and Q a 3 obtained from simulations are in good agree- 
ment with the analytical predictions #T3 = 0a3 = l/a 
[Eq. (§)]. 



0T3 



0A3 



l/a 



0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



2.29 
1.92 
1.69 
1.39 
1.26 
1.06 



2.38 
1.95 
1.71 
1.43 
1.27 
1.10 



2.50 
2.00 
1.67 
1.43 
1.25 
1.11 



Finally, our simulations show that all three crossover 
scales nix, "2 X and n3 X do not depend on the length of 
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the signal N max , since F v (n\ 
on N max as shown in Eqs. 



do not depend 
. and jis]). 



C. Higher order DFA on pure sinusoidal trend 



Since the scaling function F v (n) for correlated noise shifts 
vertically to lower values when higher order DFA-£ is used 
[see the discussion in Appendix |X] and Sec. V_B|, n 3x ex- 
hibits a slight shift to larger scales. 



In the previous Sec. IV B, we discussed how sinusoidal 



trends affect the scaling behavior of correlated noise when 
the DFA-1 method is applied. Since DFA-1 removes only 
constant trends in data, it is natural to ask how the ob- 
served scaling results will change when we apply DFA of 
order £ designed to remove polynomial trends of order 
lower than £. In this section, we first consider the rms 
fluctuation Fs for a sinusoidal signal and then we study 
the scaling and crossover properties of F^s for correlated 
noise with superimposed sinusoidal signal when higher 
order DFA is used. 

We find that the rms fluctuation function Fs does not 
depend on the length of the signal N max , and preserves 
a similar shape when different order-^ DFA method is 
used [Fig. M. In particular, Fs exhibits a crossover at a 
scale ri2x proportional to the period T of the sinusoidal: 
«2x ~ t 9t2 with 6*T2 ~ 1- The crossover scale shifts 
to larger values for higher order £ [Fig. || and Fig. ||. 
For the scale n < n,2x, Fs exhibits an apparent scaling: 



F H 



with an effective exponent as = 



1 . For 



DFA-1, we have £ = 1 and recover as — 2 as shown in 
Eq. (fjj). For n > n 2x , Fg(n) is a constant independent 
of the scale n, and of the order £ of the DFA method in 
agreement with Eq. (p"3|). 

Next, we consider F v g(n) when DFA-£ with a higher 
order I is used. We find that for all orders £, F v s(n) 
does not depend on the length of the signal N max and 
exhibits three crossovers — at small, intermediate and 
large scales — similar behavior is reported for DFA-1 in 
Fig. [| Since the crossover at small scales, ni x , and the 
crossover at large scale, 71.3 x , result from the "competi- 
tion" between the scaling of the correlated noise and the 
effect of the sinusoidal trend (Figs. [| and ^), using the 
superposition rule [Eq. fllO|)] we can estimate n\ x and 
ri3x as the intercepts of F v (n) and Fs(n) for the general 
case of DFA-1 

For nix we find the following dependence on the pe- 
riod T, amplitude ^4sj the correlation exponent a of the 
noise, and the order £ of the DFA-£ method: 



For the crossover ri2x in F v s(n) at F ri s(n) at inter- 
mediate scales, we find: ri2 X ~ T. This relation is 
independent of the order £ of the DFA and is identical to 
the relation found for Fs(n) [Eq. (fil|)]. ri2x also exhibits 
a shift to larger scales when higher order DFA is used 
[see Fig.l. 



The reported here features of the crossovers in F^s (n) 
can be used to identify low-frequency sinusoidal trends in 
noisy data, and to recognize their effects on the scaling 
properties of the data. This information may be useful 
when quantifying correlation properties in data by means 
of scaling analysis. 
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FIG. 9. Comparison of the results of different order 
DFA on a sinusoidal trend. The sinusoidal trend is given 
by the function 64sin(27ri/2 x ) and the length of the signal 
is Nmax = 2 17 . The spurious singularities (spikes) arise from 
the discrete data we use for the sinusoidal function. 



{T/M 



a/(<+i-a) 



(16) 



For DFA-1, we have 1=1 and we recover Eq. (|1J). In 
addition, nix is shifted to larger scales when higher order 
DFA-£ is applied, due to the fact that the value of Fs (n) 
decreases when £ increases (as = £ + 1, see Fig. [)]). 

For the third crossover observed in F^s (n) at large scale 
7i3x we find for all orders £ of the DFA-^ the following 
scaling relation: 



n 3 > 



(TA S ) 



l/a 



(17) 
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V. NOISE WITH POWER-LAW TRENDS 
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FIG. 10. Crossover behavior of the rms fluctuation function 
F v p(n) (circles) for correlated noise (of length N max = 2 17 ) 
with a superimposed power-law trend u(i) = Api x . The rms 
fluctuation function F v (n) for noise (solid line) and the rms 
fluctuation function Fp(n) (dash line) are also shown for com- 
parison. DFA-1 method is used, (a) F v p(n) for noise with 
correlation exponent ct\ = 0.9, and power-law trend with am- 
plitude Ap = 1000/(A r mal ) 0,4 and positive power A = 0.4; (b) 
F v p(n) for Brownian noise (integrated white noise, a\ — 1.5), 
and power-law trend with amplitude Ap = 0.01/ (N max )~ ' 7 
and negative power A = —0.7. Note, that although in both 
cases there is a "similar" crossover behavior for F v p(n), the 
results in (a) and (b) represent completely opposite situa- 
tions: while in (a) the power-law trend with positive power 
A dominates the scaling of F v p(n) at large scales, in (b) the 
power-law trend with negative power A dominates the scaling 
at small scales, with arrow we indicate in (b) a weak crossover 
in Fp(n) (dashed lines) at small scales for negative power A. 

In this section we study the effect of power-law trends 
on the scaling properties of noisy signals. We consider 



the case of correlated noise with superposed power-law 
trend u(i) = Api x , when Ap is a positive constant, 
i = 1, N max , and N max is the length of the signal. 
We find that when the DFA-1 method is used, the rms 
fluctuation function F v p (n) exhibits a crossover between 
two scaling regions [Fig. |l0| . This behavior results from 
the fact that at different scales n, either the correlated 
noise or the power-law trend is dominant, and can be 
predicted by employing the superposition rule: 



[F,p(n)]M^>)] 2 + [F P (n)f 



(18) 



where F v (ri) and Fp(n) are the rms fluctuation function 
of noise and the power-law trend respectively, and F v p (n) 
is the rms fluctuation function for the superposition of 
the noise and the power-law trend. Since the behavior of 
F^n) is known (Eq. (^|) and Appendix ^), we can un- 
derstand the features of F^pin), if we know how Fp(n) 
depends on the characteristics of the power-law trend. 
We note that the scaling behavior of F v p(n) displayed 
in Fig. |Io](a) is to some extent similar to the behavior of 
the rms fluctuation function F^^n) for correlated noise 
with a linear trend [Fig. [l| — e.g. the noise is dominant 
at small scales n, while the trend is dominant at large 
scales. However, the behavior Fp(n) is more complex 
than that of Fi,(n) for the linear trend, since the effec- 
tive exponent a\ for Fp(n) can depend on the power A 
of the power-law trend. In particular, for negative val- 
ues of A, Fp(n) can become dominated at small scales 
(Fig. |l0|(b)) while F v (ri) dominates at large scales — a 
situation completely opposite of noise with linear trend 
(Fig. [l]) or with power-law trend with positive values for 
the power A. Moreover, Fp(n) can exhibit crossover be- 
havior at small scales [Fig. |10](b)] for negative A which 
is not observed for positive A. In addition Fp(n) de- 
pends on the order £ of the DFA method and the length 
Nmax of the signal. We discuss the scaling features of 
the power-law trends in the following three subsections. 



A. Dependence of Fp(n) on the power A 



First we study how the rms fluctuation function Fp (n) 
for a power-law trend u(i) = Api x depends on the power 
A. We find that 



F P (n) - A P n° 



(19) 



where a\ is the effective exponent for the power-law 
trend. For positive A we observe no crossovers in Fp(n) 
(Fig. |l(](a)). However, for negative A there is a crossover 
in Fp(n) at small scales n (Fig. |l^(b)), and we find that 
this crossover becomes even more pronounced with de- 
creasing A or increasing the order I of the DFA method, 
and is also shifted to larger scales [Fig. |Tl](a)]. 
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FIG. 11. Scaling behavior of rms fluctuation function 
Fp(n) for power-law trends, u(i) ~ i A , where i = 1, N max 
and Nmax = 2 lr is the length of the signal, (a) For A < 0, 
Fp(n) exhibits crossover at small scales which is more pro- 
nounced with increasing the order £ of DFA-^ and decreasing 
the value of A. Such crossover is not observed for A > when 
Fp(n) ~ n° J for all scales n [see Fig. |lo|(a)]. (b) Dependence 
of the effective exponent a\ on the power A for different order 
£ — 1,2,3 of the DFA method. Three regions are observed 
depending on the order £ of the DFA: region I (A > £ — 0.5), 
where a\ m £ + 1; region II ( — 1.5 < A < £ — 0.5), where 
a\ — A + 1.5; region III (A < —1.5), where a\ « 0. We note 
that for integer values of the power A = 0, 1, ...,£ — 1, where £ 
is the order of DFA we used, there is no scaling for Fp (n) and 
ct\ is not defined, as indicated by the arrows, (c) Asymptotic 
behavior near integer values of A. Fp(n) is plotted for A — > 1 
when DFA- 2 is used. Even for A — 1 = 10~ 6 , we observe at 
large scales n a region with an effective exponent ~ 2.5, 
This region is shifted to infinitely large scales when A = 1. 

Next, we study how the effective exponent a\ for Fp(n) 
depends on the value of the power A for the power-law 
trend. We examine the scaling of Fp(n) and estimate 
ce\ for —4 < A < 4. In the cases when Fp(n) exhibits a 
crossover, in order to obtain a\ we fit the range of larger 
scales to the right of the crossover. We find that for any 
order i of the DFA-^ method there are three regions with 
different relations between a\ and A [Fig. |IT|(b)1: 

(i) a x w £ + 1 for A > £ - 0.5 (region I); 

(ii) a x w A + 1.5 for -1.5 < A < £ - 0.5 (region II); 

(iii) a\ w for A < — 1.5 (region III). 

Note, that for integer values of the power A (A = 
0, 1, — 1), i.e. polynomial trends of order m — 1, 
the DFA-£ method of order £ > m — 1 (£ is also an in- 
teger) leads to Fp(n) w 0, since DFA-^ is designed to 
remove polynomial trends. Thus for a integer values of 
the power A there is no scaling and the effective exponent 
ax is not defined if a DFA-£ method of order £ > A is used 
[Fig. [ll[] . However, it is of interest to examine the asymp- 
totic behavior of the scaling of Fp(n) when the value of 
the power A is close to an integer. In particular , we 
consider how the scaling of Fp(n) obtained from DFA-2 
method changes when A — > 1 [Fig. |TT|(c)] . Surprisingly, 
we find that even though the values of Fp(n) are very 
small at large scales, there is a scaling for Fp(n) with a 
smooth convergence of the effective exponent ax — ► 2.5 
when A — > 1, according to the dependence ax ~ A + 1.5 
established for region II [Fig. [ll](b)]. At smaller scales 
there is a flat region which is due to the fact that the 
fluctuation function Y(i) (Eq. (||)) is smaller than the 
precision of the numerical simulation. 



B. Dependence of Fp(n) on the order £ of DFA 

Another factor that affects the rms fluctuation func- 
tion of the power-law trend Fp(n), is the order I of the 



DFA method used. We first take into account that: 



(c) 



(1) for integer values of the power A, the power-law 
trend u{i) — Api x is a polynomial trend which 
can be perfectly filtered out by the DF A met hod 
of o rder t > A, and as discussed in Sec. Ill B| and 
Sec. VA [see Fig.|ll](b) and (c)], there is no scaling 
for Fp(n). Therefore, in this section we consider 
only non-integer values of A. 



(2) for a given value of the power A, the effective ex- 
ponent a\ can take different values depending on 
the order £ of the DFA method we use [see Fig. Q 
- e.g. for fixed A > I - 0.5, a x ^ £ + 1. There- 
fore, in this section, we consider only the case when 
X < £-0.5 (Region II and III). 
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FIG. 12. Effect of higher order DFA-^ on the rms fluctua- 
tion function F v p(n) for correlated noise with superimposed 
power-law trend, (a) F v p(n) for anticorrelated noise with 
correlation exponent a = 0.1 and a power-law u(i) — Api x , 
where A P = 25/ (N max ) 0A , N max = 2 17 and A = 0.4. Results 
for different order I = 1,2, 3 of the DFA method show (i) a 
clear crossover from a region at small scales where the noise 
dominates F v p(n) ~ F v (n), to a region at larger scales where 
the power- law trend dominates F v p(n) « Fp(n), and (ii) a 
vertical shift A in F^p with increasing I. (b) Dependence 
of the vertical shift A in the rms fluctuation function Jp(n) 
for power-law trend on the order I of DFA-^ for different val- 
ues of A: A ~ £ tI ~ x \ We define the vertical shift A as the 
y- intercept of Fp(n): A = Fp(n = 1). Note, that we consider 
only non-integer values for A and that we consider the region 
A < I — 0.5. Thus, for all values of A the minimal order i 
that can be used in the DFA method is I > A + 0.5. e.g. for 
A = 1.6 the minimal order of the DFA that can be used is 
I = 3 (for details see Fig. |lT|(b)), (c) Dependence of r on the 
power A (error bars indicate the regression error for the fits of 
A(Z) in (b)). (d) Comparison of t(qa) for Fp(n) and t(q) for 
F v (n). Faster decay of r(a\) indicates larger vertical shifts 
for Fp(n) compared to F v (n) with increasing order I of the 
DFA-1 
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Since higher order DFA-^ provides a better fit for the 
data, the fluctuation function Y(i) (Eq. (|j)) decreases 
with increasing order £. This leads to a vertical shift 
to smaller values of the rms fluctuation function F(n) 
(Eq. ([!])). Such a vertical shift is observed for the rms 
fluctuation function F v (n) for correlated noise (see Ap- 
pendix as well as for the rms fluctuation function of 
power-law trend Fp(n). Here we ask how this vertical 
shift in F v (n) and Fp(n) depends on the order £ of the 
DFA method, and if this shift has different properties 
for F v (n) compared to Fp (n) . This information can help 
identify power-law trends in noisy data, and can be used 
to differentiate crossovers separating scaling regions with 
different types of correlations, and crossovers which are 
due to effects of power-law trends. 

We consider correlated noise with a superposed power- 
law trend, where the crossover in F v p (n) at large scales n 
results from the dominant effect of the power-law trend 
- F ri p(n) « F P (n) (Eq. @ and Fig. [Ufa)). We choose 
the power A < 0.5, a range where for all orders i of the 
DFA method the effective exponent ct\ of Fp (n) remains 
the same — i.e. a\ = A +1.5 (region II in Fig. [il|(b)). For 
a superposition of an anticorrelated noise and power-law 
trend with A = 0.4, we observe a crossover in the scaling 
behavior of F r) p(n), from a scaling region characterized 
by the correlation exponent a = 0.1 of the noise, where 
i 7 ' I) p(n) w F r/ (n), to a region characterized by an effective 
exponent oi\ — 1.9, where F rj p(n) Fp (n) , for all orders 
t = 1,2,3 of the DFA-£ method [Fig. |lf(a)]. We also 
find that the crossover of F v p(n) shifts to larger scales 
when the order I of DFA-£ increases, and that there is 
a vertical shift of F v p(n) to lower values. This vertical 
shift in F v p(n) at large scales, where F v p(n) = Fp(n), 
appears to be different in magnitude when different or- 
der I of the DFA-£ method is used [Fig. ||(a)]. We also 
observe a less pronounced vertical shift at small scales 
where F^pin) « F r/ (n). 

Next, we ask how these vertical shifts depend on the 
order i of DFA-£. We define the vertical shift A as the 
y-intercept of Fp(n): A = Fp(n = 1). We find that the 
vertical shift A in Fp(n) for power-law trend follows a 
power law: A ~ £ T ^ X \ We tested this relation for orders 
up to I — 10, and we find that it holds for different val- 
ues of the power A of the power-law trend [Fig. |l2](b)]. 
Using Eq. ( |l9| ) we can write: Fp(n)/Fp(n = 1) = n ax , 
i.e. Fp(n) ~ F P {n = 1). Since F P (n = 1) = A ~ £ T ^ 
[Fig. P(b)], we find that: 



F P (n) ~fW. 



(20) 



We also find that the exponent r is negative and is a 
decreasing function of the power A [Fig. [12|(c)]. Because 
the effective exponent a a which characterizes Fp(n) de- 
pends on the power A [see Fig. |H|(b)], we can express the 
exponent r as a function of a\ as we show in Fig. [l^(d) . 
This representation can help us compare the behavior of 
the vertical shift A in Fp(n) with the shift in F v (n). For 
correlated noise with different correlation exponent a, we 



observe a similar power-law relation between the vertical 
shift in F v (n) and the order £ of DFA-£: A ~ £ T{ - a \ where 
t is also a negative exponent which decreases with a. In 
Fig. |l2|(d) we compare r(a\) for Fp(n) with r(a) for 
F^in), and find that for any a\ — a, r(a\) < r(a). This 
difference between the vertical shift for correlated noise 
and for a power-law trend can be utilized to recognize 
effects of power-law trends on the scaling properties of 
data. 



C. Dependence of Fp(n) on the signal length N max 



Here, we study how the rms fluctuation function _Fp (n) 
depends on the length N max of the power-law signal 
u(i) = Api x (i = 1, N max ). We find that there is a 
vertical shift in Fp (n) with increasing N max [Fig. |l3|(a)] . 
We observe that when doubling the length N max of the 
signal the vertical shift in Fp(n), which we define as 
Fp Nma * /Fornax ^ remams t ne same, independent of the 
value of N max . This suggests a power-law dependence of 
Fp(n) on the length of the signal: 



F P (n) ~ {N max )\ (21) 



where 7 is an effective scaling exponent. 

Next, we ask if the vertical shift depends on the power 
A of the power-law trend. When doubling the length 
N m ax of the signal, we find that for A < £ ~ 0.5, where £ 
is the order of the DFA method, the vertical shift is a con- 
stant independent of A [Fig. [l^(b)]. Since the value of the 
vertical shift when doubling the length N max is 2 7 (from 
Eq. (21)), the results in Fig. [l^(b) show that 7 is inde- 
pendent of A when A < £— 0.5, and that — log2 7 w —0.15, 
i.e. the effective exponent 7 « —0.5. 

For A > £ — 0.5, when doubling the length N max of 
the signal, we find that the vertical shift 2 7 exhibits the 
following dependence on A: — log 10 2 7 = log 10 2 A ~ £ , and 
thus the effective exponent 7 depends on A — 7 = A — £. 
For positive integer values of A (A = £), we find that 
7 = 0, and there is no shift in Fp(n), suggesting that 
Fp(n) does not depend on the length N max of the signal, 
when DFA of order £ is used [Fig. [Dj| . Finally, we note 
that depending on the effective exponent 7, i.e. on the 
order £ of the DFA method and the value of the power A, 
the vertical shift in the rms fluctuation function Fp(n) 
for power-law trend can be positive (A > £), negative 
(A < £), or zero (A = £). 
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and the effect of the length of the signal N max (Eq. (pl|)), 
we reach the following expression for the rms fluctuation 
function Fp(n) for a power-law trend u(i) = Api x : 



Fp(n)~A P -n a * ■rW-(7V max ) 



7(A) 



(22) 



For correlated noise, the rms fluctuation function F v (n) 
depends on the box size n (Eq. (^)) and on the order I 
of DFA-^ (Sec. [VB| and Fig. |§(a), (d)), and does not 
depend on the length of the signal N max . Thus we have 
the following expression for 



(23) 



To estimate the crossover scale n x observed in the 
apparent scaling of F v p(n) for a correlated noise su- 
perposed with a power-law trend [Fig. 0(a), (b) and 
Fig. [I^(a)], we employ the superposition rule (Eq. (|l8|)). 
From Eq. ( p2|) and Eq. fl25|), we obtain n x as the inter- 
cept between Fp(n) and F^in): 



Al r(X)-T(a) {Nmax) - 



1/ {a— a\) 



(24) 



To test the validity of this result, we consider the case of 
correlated noise with a linear trend. For the case of a lin- 
ear trend (A = 1) when DFA-1 (£ = 1) is applied, we have 
ax = 2 (see Appendix |c| and Sec. |VA| , Fig.]l|(b)). Since 
in this c ase A = £ — 1 > £ — 0.5 we have 7 — A — £ = 
(see Sec.VC Fig. |l3|(b)), and from Eq. (|j) we recover 
Eq. @. 



FIG. 13. Dependence of the rms fluctuation function Fp(n) 
for power-law trend u(i) = Api x , where i = 1, N max , on 
the length of the trend N max . (a) A vertical shift is ob- 
served in Fp(n) for different values of N ma x — Ni max and 
N2max- The figure shows that the vertical shift , defined as 



F p lmax (n)/Fp 2max (n), does not depend on N max but only on 
the ratio N lmax /N 2m ax, suggesting that Jp(n) ~ (jV maa; ) 7 . 
(b) Dependence of the vertical shift on the power A. For 
A < i — 0.5 (I is the order of DFA), we find a flat (constant) 
region characterized with effective exponent 7 = —0.5 and 
negative vertical shift. For A > ^ — 0.5, we find an exponential 
dependence of the vertical shift on A. In this region, 7 = A — £, 
and the vertical shift can be negative (if A < £) or positive (if 
A > 1). the slope of -log 10 (Fp Nma * {n) / Fp ma * {n)) vs. A is 
— log 10 2 due to doubling the length of the signal Nmax- This 
slope changes to — log 10 m when N max is increased m times 
while 7 remains independent of Nmax- For A = t there is 
no vertical shift, as marked with x. Arrows indicate integer 
values of A < £, for which values the DFA-£ method filters out 
completely the power-law trend and _Fp = 0. 



D. Combined effect on Fp(n) of A, I and N max 

We have seen that, taking into account the effects of 
the power A (Eq. ©), the order I of DFA-£ (Eq. (§C 



VI. CONCLUSION AND SUMMARY 

In this paper we show that the DFA method performs 
better than the standard R/S analysis to quantify the 
scaling behavior of noisy signals for a wide range of cor- 
relations, and we estimate the range of scales where the 
performance of the DFA method is optimal. We consider 
different types of trends superposed on correlated noise, 
and study how these trends affect the scaling behavior of 
the noise. We demonstrate that there is a competition be- 
tween a trend and a noise, and that this competition can 
lead to crossovers in the scaling. We investigate the fea- 
tures of these crossovers, their dependence on the proper- 
ties of the noise and the superposed trend. Surprisingly, 
we find that crossovers which are a result of trends can 
exhibit power-law dependences on the parameters of the 
trends. We show that these crossover phenomena can be 
explained by the superposition of the separate results of 
the DFA method on the noise and on the trend, assum- 
ing that the noise and the trend are not correlated, and 
that the scaling properties of the noise and the appar- 
ent scaling behavior of the trend are known. Our work 
may provide some help to differentiate between differ- 
ent types of crossovers — e.g. crossovers which separate 
scaling regions with different correlation properties may 
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differ from crossovers which are an artifact of trends. The 
results we present here could be useful for identifying the 
presence of trends and to accurately interpret correlation 
properties of noisy data. 
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APPENDIX A: NOISE 



The standard signals we generate in our study are un- 
corrected, correlated, and anticorrelated noise. First we 
must have a clear idea of the scaling behaviors of these 
standard signals before we use them to study the effects 
from other aspects. We generate noises by using a mod- 
ified Fourier filtering method |6j|. This method can ef- 
ficiently generate noise, u(i) (i — 1, 2, 3, A max ), with 
the desired power-law correlation function which asymp- 

i+t 

totically behaves as: < | ^ M (j)| 2 >^ t 2a ■ By default, a 

generated noise has standard deviation a = 1. Then we 
can test DFA and R/S by applying it on generated noises 
since we know the expected scaling exponent a. 
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FIG. 14. Scaling behavior of noise with the scaling ex- 



ponent a. The length of noise N n 



= 2 1 



(a) Rescaled 



range analysis (R/S) (b) Order 1 detrended fluctuation anal- 
ysis (DFA-1) (c) Order 2 detrended fluctuation analysis. We 
do the linear fitting for R/S analysis and DFA-1 in three re- 
gions as shown and get ai, 0:2 and 0:3 for estimated a, which 
are listed in the Table|vJ and Table.0. We find that the 
estimation of a is different in the different region. 

Before doing that, we want to briefly review the algo- 
rithm of R/S analysis. For a signal u(i)(i = 1, A max ), 
it is divided into boxes of equal size n. In each box, 
the cumulative departure, Xi (for fc-th box, i = kn + 
1, kn + n), is calculated 



X i = £ M-/)- < U >) 
j—kn-\-l 



(Al) 



(k+l)n 

where < u >= ^2 U W > an( i the rescaled range 

i—kn+l 

R/S is defined by 
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R/S = S- 



where S 



max Xi 

kn+l<i<(k+l)n 



min Xi 

kn+l<i<{k+l)n 



study shows that the scaling exponent a depends on scale 
n. The estimated a is different for the different regions 



(A2) 



( u (j)~ < u >) 2 i s the standard de- 

viation in each box. The average of rescaled range in 
all the boxes of equal size n, is obtained and denoted by 
< R/S >. Repeat the above computation over different 
box size n to provide a relationship between < R/S > 
and n. According to Hurst's experimental study J64J , a 
power-law relation between < R/S > and the box size n 
indicates the presence of scaling: < R/ S >~ n a . 

Figure [u] shows the results of R/S, DFA-1 and DFA- 
2 on the same generated noises. Loosely speaking, we 
can see that F(n) (for DFA) and R/S (for R/S analysis) 
show power-law relation with n as expected: F(n) ~ n a 
and R/S ~ n a . In addition, there is no significant dif- 
ference between the results of different order DFA except 
for some vertical shift of the curves and the little bend- 
down for small box size n. The bent-down for very small 
box of F(n) from higher order DFA is because there are 
more variables to fit those few points. 

TABLE IV. Estimated a of correlation noise from 
R/S analysis in three regions as shown in Figjl4|(a). a 
is the input value of the scaling exponent, ai is the 
estimated in the region 1 for 4 < n < 32, 012 in the 
region 2 for 32 < n < 3162 and 03 in the region 3 for 
3126 < n < 2 17 . Noise are the same as used in Table.j^. 



a 


ai 


a 2 


«3 


0.1 


0.44 


0.23 


0.12 


0.3 


0.52 


0.37 


0.23 


0.5 


0.62 


0.52 


0.47 


0.7 


0.72 


0.70 


0.45 


0.9 


0.81 


0.87 


0.63 



TABLE V. Estimated a of correlation noise from 
DFA-1 in three regions as shown in Figjl4|(b). a is 
the input value of the scaling exponent, al is the es- 
timated in the region 1 for 4 < n < 32, a2 in the 
region 2 for 32 < n < 3162 and a3 in the region 3 for 
3126 < n < 2 17 . 



a 


Otl 


Q2 


az 


0.1 


0.28 


0.15 


0.08 


0.3 


0.40 


0.31 


0.22 


0.5 


0.55 


0.50 


0.35 


0.7 


0.72 


0.69 


0.55 


0.9 


0.91 


0.91 


0.69 



Ideally, when analyzing a standard noise, F(n) (DFA) 
and R/S (R/S analysis) will be a power-law function 
with a given power: a, no matter which region of F(n) 
and R/S is chosen for calculation. However, a careful 



of F(n) and R/S as illustrated by Figs. 14(a) and |l4|(b) 
and by Tables IV and It is very important to know the 
best fitting region of DFA and R/S analysis in the study 
of real signals. Otherwise, the wrong a will be obtained 
if an inappropriate region is selected. 



In order to find the best region, we first determine the 
dependence of the locally estimated a, ot\ oc , on the scale 
n. First, generate a standard noise with given scaling 
exponent a; then calculate F(n) (or R/S), and obtain 
&\oc(n) by local fitting of F(n) (or R/S). Same random 
simulation is repeated 50 times for both DFA and R/S 
analysis. The resultant average ai oc (n), respectively, are 
illustrated in Figjl^ for DFA-1 and R/S analysis. 



If a scaling analysis method is working properly, then 
the result a\ oc (n) from simulation with a would be a 
horizontal line with slight fluctuation centered about 
ct\oc(n) = a. Note from Fig.[l5] that such a horizontal 
behavior does not hold for all the scales n but for a cer- 
tain range from n m j n to n max . In addition, at small scale, 
R/S analysis gives aq c > a if a < 0.7 and a\ oc < a if 
a > 0.7, which has been pointed out by Mandelbrot pq ] 
while DFA gives a\ oc > a if a < 1.0 and a\ oc < a if 
a > 1.0. 



It is clear that the smaller the n m i n and the larger the 
"max, the better the method. We also perceive that the 
expected horizontal behavior stops because the fluctua- 
tions become larger due to the under-sampling of F(n) 
or R/S when n gets closer to the length of the signal 
A max . Furthermore, it can be seen from Fig.|l5| that 
"max ~ jg -^max independent of a (if the best fit region 
exists), which is why one tenth of the signal length is the 
maximum box size when using DFA or R/S analysis. 



1.1 ™, — — — — | 

(a) R/SN =2 14 

max . ••••• • 

0.9 . • ••• • 

+ + + + + + + + + + + + + + + + + + + + T3 

+ - £ 
ra 

0.7 - «-Xxrx-xi*xx-xXiK->:- !< - < -< !( <D 

O T3 

* * * * - -2 

o ****** 5 
a" 0.5 oa=0.1 <rv *_***_*_*_****-«--^£ 

3-* a = - 5 ^o,__^o^_<>^ - 3 
U ' J < a=0.7 o o n ! § 

-+a=0.9 °?°°°oo 00 - 
0.1 • ink\«nitcil </.-(). I { 

10° 10 1 10 2 n min io 3 10 4 
n 



16 



1.1 



0.9 



0.7 

u 

d"o.5 

0.3 
0.1 



FIG. 15. The estimated a from local fit (a) R/S anal- 



0.1 



(b) R/S N ma -2 2 :...............«'... - 

+ ■ ® 

********* 

a=0.1 o d ****|***jt**************---^, 2 

a=0.3*°°^ o : . § 

«= - 5 _ -^°o CT -- ^^^c^o.--- 3 
01=0.7 °°o 0o ! o 'f 

a = - 9 30 ? 00000o °ooooooooo 0o " " 
integrated <3l=0T ; ---- 

10° 10 2 n min 10 4 

n 



(c) DFA-1 N =2 

v ' max 



max 



1.1 

0.9 

0.7 - X--x--x-x-x-x--x-x-x-x--x-X-*-x-^ 

i 

0.5 - o oc=0.1 ***♦♦»♦»»♦*♦»♦< 



a=0.3 



* * * 
* * 



0.3 



a=0.5 _°__ 0o *_o_<>o*_< 



x a=0.7 ° o , 
+ a=0.9 
• integrated a=6.1 



o o 



Oo poooooo 0< 



10" 



10 1 n min10 2 

n 



10 J 



10' 



1.1 

0.9 

0.7 

) 

; 

0.5 

0.3 

0.1 
10 



(d) DFA-1 N =2 



20 



+-H-HH-+++++++ ,, + T ------ 



^xx^*sxxxm*^xx^-x« x ^-< ^- - 



a=0.3 
a=0.5 



a=0.9 10 °°?oc 
• integrated a=0.1 ; 

1 10 1 n 



°o«X)oqooo! 



CfXxQOQQQQQOO- - - 



"O 

w 
ca 

CD 

O CD 



10 J 



10= 



n 



ysis, the length of signal JVn 



(b)R/S analysis, 



iVmax = 2 ZU . (c) DFA-1, AW = 2 14 (d) DFA-1, AW = 2 M . 
ot\ oc come from the average of 50 simulations. If a technique 
is working, then the data for scaling exponent a should be a 
weakly fluctuating horizontal line centered about a\ oc = a. 
Note that such a horizontal behavior does not hold for all the 
scales. Generally, such a expected behavior begins from some 
scale n m ; n , holds for a range and ends at a larger scale n max . 
For DFA-1, n m ; n is quite small a > 0.5. For R/S analysis, 
^min is small only when a ~ 0.7. 
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FIG. 16. The starting point of good fit region, ra m i n , for 
DFA-1 and R/S analysis. The results are obtained from 50 
simulations, in which the length of noise is AW = 2 20 . The 
condition for a good fit is Aa = |qi oc — a \ < 0.01. The data 
for a > 1.0 shown in the shading area are obtained by apply- 
ing analysis on the integrations of noises with a < 1.0. It is 
clear that DFA-1 works better than R/S analysis because its 
n min is always smaller than that of R/S analysis. 

On the contrary, n m i n does not depend on the iV max 
since aiocl") at small n hardly changes as N max varies 
but it does depend on a. Thus, we obtain n min quan- 
titatively as shown in Fig.|l^. For R/S analysis, only 
for a ~ 0.7, n m ; n is small; for a a little away from 0.7 
(for example, 0.5), n m i n becomes very large and close to 
"■max, indicating that the best fit region will vanish and 
R/S analysis does not work at all. Comparing to R/S, 
DFA works better since n m i n is quite small for a > 0.5 
correlated signals. 

One problem remains for DFA, n m j n for small a (< 0.5) 
is still too large comparing to those for large a(> 0.5). 
We can improve it by applying DFA on the integration 
of the noise with a < 0.5. The resultant new expected 
a for the integrated signal would be a Q = a + 1, while 
the n m ; n for the integrated signal becomes much smaller 
as shown also in Fig.[l6](shading area a > 1). Therefore, 
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for a noise with a < 0.5, it is best to estimate the scaling 
exponent a of the integrated signal first and then obtain 
a by a — a — 1. This is what we did in the following 
sections to those anticorrelated signals. 



APPENDIX B: SUPERPOSITION LAW FOR DFA 

For two uncorrelated signals f(i) and g(i), their root 
mean square fluctuation functions are Ff(n) and F g (n) 
respectively. We want to prove that for the signal 
f(i) + g(i), its rms fluctuation 



F 



' f+g (n) = ^F f (nr + F g (nf 



(Bl) 



Consider three signals in the same box first. The in- 
tegrated signals for /, g and f + g are ?//(«), Ug{i) and 
yf+ g (i) and their corresponding trends are y^ lt , y g rt ,y^ g 
(i = 1,2, ...,n, n is the box size). Since yf +g (i) — 
yf(i)+y g (i) and combine the definition of detrended fluc- 
tuation function Eq.^|, we have that for all boxes 



Y f+g (i) = Y f (i)+Y g {i), 



(B2) 



where Yf +g is the detrended fluctuation function for the 
signal f + g, Yf(i) is for the signal / and Y g (i) for g. Fur- 
thermore, according to the definition of rms fluctuation, 
we can obtain 



*>+,(»> 



\ 



max , 
i—l 



Nmax 



\ 



N„ 



N„ 



£ [Y f (i)+Y g (i)f 



i=i 



where £ is the number of boxes and k means the kth box. 
If / and g are not correlated, neither are Yf(i) and Y g (i) 
and. thus, 



E 



Y f (i)Y g (i) = 0. 



(B4) 



From Eq.B4 and Eq.|B4|, we have 

*>+,(») 



\ 



1 



N„ 



[Ff{r 



(B5) 



APPENDIX C: DFA-1 ON LINEAR TREND 

Let us suppose a linear time series u(i) = A^i. The 
integrated signal is 



3 = 1 



(CI) 



Let as call N max the size of the series and n the size of 
(B3) the box. The rms fluctuation F^(n) as a function of n 

and Nm.n-r. is 



F L (n) = A L 



\ 



N mam /n kn 

E E 

fc=l i=(fe-l)n+l 



(a k + b k i) 



(2) 



where a k and b k are the parameters of a least-squares 
fit of the fc-th box of size n. a k and b k can be determined 
analytically, thus giving: 



(ii. = 1 — n 2 + -n 2 k + —n — -k 2 n 2 

12 2 12 2 



(3) 



1 



1 



bk = 1 — —n + kn + — 



2 2 

With these values, F^(n) can be evaluated analytically 
1 



(4) 



F L (n) = A L — V(5n 4 + 25n 3 + 25n 2 - 25n - 30) (5) 

The dominating term inside the square root is 5n and 
then one obtains 



(6) 



leading directly to an exponent of 2 in the DFA. An im- 
portant consequence is that, as F(n) does not depend on 



Nmax, f° r linear trends with the same slope, the DFA 
must give exactly the same results for series of different 
sizes. This is not true for other trends, where the expo- 
nent is 2, but the factor multiplying n can depend on 

N max - 



APPENDIX D: DFA-1 ON QUADRATIC TREND 

Let us suppose now a series of the type u(i) = Aq{ 2 . 
The integrated time series y(i) is 



2i 6 + 3i 2 + i 
6 



(1) 



As before, let us call N m ax and n the sizes of the se- 
ries and box, respectively. The rms fluctuation function 
Fq(ji) measuring the rms fluctuation is now defined as 
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Fq(n) = Ac, 



\ 



N, 



. N maI /n kn 

— E E 

nax k=l i=(fe-l)n+l 



2i 3 + 3i 2 + i 
6 



- (a* + M) 



(2) 



where and 6fc are the parameters of a least-squares fit of the k-th box of size n. As before, a k and 6^ can be 
determined analytically, thus giving: 



a k = — n 3 + n 3 k 2 - — n 3 fc + — n 2 k - —n 2 + —n - -k 3 n 3 - -n 2 k 2 + —hi 
15 15 30 60 20 3 2 15 



(3) 



i J 2 

6 fe = 



21.2 



n k 



2, , 2 1 
n k + kn n H — - 



Once afc and are known, F(n) can be evaluated, giving: 



10 



(4) 



F Q (n) = Aq^V-21 {n A + 5n 3 + 5n 2 - 5ra - 6) (32? 



, 2 - 6n - 81 - 210iV„ 



140iV 2 



(5) 



As N max > n, the dominant term inside the square 
root is given by U0N^ ax x 21n 4 = A Q 2940n 4 N^ ax , and 
then one has approximately 



F Q (n) n A Q — ^/2940nM- 



2 

max 



Ac 



'90 



15V„ 



(6) 



leading directly to an exponent 2 in the DFA analysis. 
An interesting consequence derived from Eq. (^) is that, 
Fq(u) depends on the length of signal N max , and the 
DFA line (logi*Q(n) versus logn) for quadratic series 
u(i) — Aqi 2 of different N max DO NOT overlap (as it 
happened for linear trends). 
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